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Запропоновано математичну модель, яка дозволяє враховувати теплову інерцію земель 
сільськогосподарського призначення в інфрачервоному діапазоні. Проаналізовано можливі шляхи 
застосування тепловізорів при контролі земель сільськогосподарського призначення. Проведено 
порівняльний аналіз даних, отриманих за допомогою математичної моделі, з даними натурного 
експерименту. 
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Вступ 

У сучасному світі все більшого значення набуває аграрний сектор економіки. 
Найбільші економіки світу мають обмежені ресурси щодо вирощування аграрних 
культур і виробництва продуктів харчування. Так, у США, більшість земель, на яких 
можливе вирощування аграрних культур, вже задіяні у цій сфері, а значне їх 
збільшення неможливе у зв'язку із об'єктивними причинами (посушливі райони 
місцевості, гори, болота і т.п.). У Китаї та Індії, де фактично проживає половина 
населення Землі, відчувається значна нестача земель, придатних для вирощування 
власної аграрної продукції. Вочевидь, аграрний сектор економіки поступово набуває 
все більшого значення. 

Особливе місце на сучасному етапі розвитку аграрного сектору займає питання 
моніторингу земель сільськогосподарського призначення (ЗСП). Це дозволяє вчасно 
реагувати на їх стан, а відтак, 1 прогнозувати подальші зміни врожайності. Одним із 
найбільш дешевих і ефективних шляхів видового моніторингу ЗСП є використання 
безпілотних літальних апаратів (БЛА). Такі БЛА обладнуються телевізійними камерами 
спостереження, що працюють у видимому та інфрачервоному (ІЧ) діапазоні електро- 
магнітного спектру. Саме ІЧ діапазон має ряд переваг, що дозволяють більш ефективно 
проводити моніторинг ЗСП. 

Таким чином, метою статті є прогнозування циклічних добових змін температури 
ЗСП для підвищення достовірності результатів моніторингу під час застосування БЛА, 
обладнаних тепловізорами. 

Постановка проблеми 

Одним із можливих способів збільшення об'ємів вирощування аграрних культур є 
інтенсивний. Він дозволяє збільшувати врожайність аграрних культур при незмінній площі 
ЗСП. Проте, для отримання максимальної врожайності, необхідно своєчасно проводити 
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моніторинг ЗСП. У даному контексті найбільш доцільно використовувати системи 
моніторингу із елементами штучного інтелекту. Це дозволить вчасно виявляти стан 
зрілості аграрних культур, проводити обробку від шкідників та хвороб і збирати врожай. 

На даний час існує ряд систем моніторингу ЗСП із елементами штучного 
інтелекту, що без контролю оператора здатні проводити елементи посівної компанії, 
здійснювати контроль за ростом аграрних культур та збирати врожай. У дані системи 
входять: сервери для обробки інформації з моніторингу ЗСП (програмні алгоритми із 
елементами штучного інтелекту), БЛА для видового контролю, автоматизована 
сільськогосподарська техніка для посіву та збору врожаю. 

Аналіз останніх досліджень і публікацій (1-4| вказує на перспективність 
досліджень у напрямку визначення циклічних добових змін температури ЗСП. Так, 
робота (1| присвячена результатам експериментів з визначення теплового контрасту 
людини на фоні грунту (рослинності) при різній фоновій обстановці протягом доби в 
діапазонах 1,9 - 5,5 17,6- 11,8 мкм. 

Існуюча математична модель опису циклічної добової зміни температури трунтів 
наведена в патенті |2|, також у даній роботі проводиться прогнозування значень 
температури протягом доби та інтервалів часу доби, коли циклічна добова зміна 
теплового контрасту максимальна та мінімальна. Наведена в роботі |2| математична 
модель не враховує всі фактори, що впливають на циклічну добову зміну температури 
ЗСП, а це, як наслідок, не дозволяє в повній мірі виконувати поставлені завдання. Перш 
за все, це викликано тим, що обмеження, які бралися при створенні математичної 
моделі, не враховували інерційних властивостей різних видів грунтів (зміну амплітуди 1 
фази коливання температури). По друге, екстремуми значень температури досягаються 
в 0:00 та о 12:00, що не в повній мірі відповідає дійсності. А вплив географічного 
фактору і радіаційного випромінювання Сонця розкриті недостатньо. Основною 
причиною цього є неврахування математичною моделлю дійсного часу досягнення 
максимального (мінімального) значення температури навколишнього середовища, 
відповідно до прогнозу погоди. У працях |3, 4| викладено основи і визначення 
циклічних добових змін теплообмінних процесів між грунтом і навколишнім 
середовищем. 


Аналіз робіт (1-4 свідчить про те, що питання прогнозування циклічних добових 
змін температури ЗСП розкрите не в повному обсязі та має ряд недоліків, тобто 
потребує додаткового і більш конкретного вивчення. 


Мета дослідження 

Таким чином, метою статті є прогнозування циклічних добових змін температури 
ЗСП для підвищення достовірності результатів моніторингу під час застосування БЛА 
обладнаних тепловізорами. 

Виклад основного матеріалу 

Перевагами використання тепловізорів під час видового моніторингу ЗСІ є 
всепогодність їх роботи (за наявності густого туману) та можливість цілодобового 
моніторингу (вдень і вночі). Дані переваги викликані особливостями розповсюдження 
електромагнітних хвиль в ІЧ діапазоні. Проте, відсутність розуміння фізичних процесів 
у навколишньому середовищі не дозволяє повною мірою використовувати наявні 
переваги тепловізорів. Протягом доби значення температури ЗСП зазнають циклічних 
змін. Це призводить до того, що в одні інтервали часу доби використання тепловізорів 
раціональне (за температурним розрізненням) (рис. 1), а в інші інтервали часу 
використання тепловізорів нераціональне (температура ЗСП та грунтових автодоріг 
практично однакові). 
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Рис. 1. Тепловізійне зображення ЗСП за допомогою БЛА 


Вище наведені особливості можна врахувати, використовуючи загальноприйняті 
методи моделювання теплообмінних процесів. Зміна температури навколишнього 
середовища та інтенсивність радіаційного випромінювання Сонця є циклічними 
процесами, тому, відповідно до теорії теплообміну, це задача нестаціонарної 
теплопровідності з граничними умовами третього роду. Існуючі наукові роботи 
дозволяють визначати інерційні властивості ЗСП, а також - кількість тепла, що 
акумулюється в них протягом півперіоду коливання температури (4-6). Проте, для 
спрощення математичних розрахунків, задається зміна температури навколишнього 
середовища у вигляді синусоїдальної (косинусоїдальної) функції, географічні фактори 
та вплив радіаційного випромінювання Сонця не враховуються. 

Знаходження розв'язку диференційного рівняння теплопровідності, з врахуванням 
всіх необхідних обмежень, є досить затратним у часі та не завжди має аналітичний 
розв'язок. Існує необхідність спрощення поставленої задачі. Це можливо за умови 
знаходження виключно зміни температури поверхні ЗСП, при циклічному добовому 
впливі на них навколишнього середовища і радіаційного випромінювання Сонця. При 
цьому, розв'язання задачі на визначення прогріву ЗСП в глибину та акумуляцію ними 
тепла не беруться до уваги (рис. 2). 


зби 
Конвекційний теплообмін 94 
У | Власне радіаційне 

11) випромінювання ЗСП. 


Радіаційне випромінювання Сонця 


Рис. 2. Умовний поділ ЗСП на частини та влив навколишнього середовища на них 
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Задача на зміну температури поверхні ЗСП зводиться до відомої задачі із 
визначення нагрівання об'єкта із нескінченно великою теплопровідністю. Характерною 
особливістю такої задачі є те, що об'єкт вважається високотеплопровідним (2-00, 
Ві-»0), тобто температура по всьому об'єму в кожний момент часу буде практично 


однаковою. У такому випадку УТ -0, а диференційне рівняння теплопровідності 
(Фур'є-Кірхгофа) набуває наступного вигляду (51: 
ФО В 


а ер, (1 

де Т - температура поверхні ЗСП, К; / - час перебування ЗСП під впливом 
теплообміну, с; Р, - теплонапруженість, Вт/м); с - питома теплоємність ЗСП, Дж/(кг-К); 
р - густина ЗСП, кг/м.. 

Негативним у рівнянні (1) є його залежність від об'єму. Для великих значень 
об'ємів температура ЗСП буде практично незмінною, тільки зі зменшенням об'єму 
вплив ставатиме більшим. Отже, значення об'єму для знаходження циклічної зміни 
температури поверхні ЗСП має бути максимально малим. 

Врахуємо одночасний вплив конвекційного і радіаційного теплообміну з 
навколишнім середовищем наступним чином: 


1 і 
щ. ех б М б Р |аб(Т, -1 Т) ра «а 510 Ук, (віп ІД б Ки 50 У» ) 
й 10) 
де Ре, Р - теплонапруженість під час конвекційного і радіаційного теплообміну, 
Вт/м3; У - об'єм грунту, м3; а - коефіцієнт теплообміну, Вт/(м?:К); 5 - площа поверхні 
ЗСП, яка бере участь у конвекційному і радіаційному теплообміні, м?; Т, - температура 
навколишнього середовища, К; а(2.) - спектральний коефіцієнт поглинання теплового 
випромінювання; 4. - сонячна постійна, (4 - 1366) Вт/м?) Тати - коефіцієнт пропускання 
атмосфери |7|; У, - умова наявності Сонця; у - кут максимальної висоти Сонця над 
горизонтом для визначеної географічної широти, град |З|; у; - кут, що враховує 
обертання Землі навколо своєї осі, град; у» - кут, що враховує обертання Землі навколо 
Сонця, град; Кк; - поправочний коефіцієнт для обертання Землі навколо своєї осі; ХК, - 
поправочний коефіцієнт для географічної широти. 

Рівняння (2) не дає змоги визначати більшість його складових, тому існує 
необхідність у їхньому поелементному описі. Так, температура середовища Т, 
описується рівнянням |21: 

Т.З То, У ОАТ, іпр' Сг), 0) 


де ДТ, - різниця між максимальним і мінімальним значеннями температури 
навколишнього середовища, К; Т,, - середня температура навколишнього середовища 
протягом доби, К. 


ре) --х Ло, 
р НН Р Р | Са 


й : Я - І І 
У рівнянні (3) введені такі позначки: ел і зай 


де й - час початку теплообміну, с; І, - час сходу Сонця (г, - 9:3,6-107), с; Ід - 
довжина світлової частини доби (/,д- 12:3,6:103), с. 
Час кореляції /хх ВИЗНАачається за допомогою рівняння 


Їк з 3,6 | 107 СЛУ З 12) (4) 


Де Ісктах - Час дОСЯГНення максимального значення температури середовища 
відповідно до прогнозу погоди, с. 
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Для знаходження у використовується наступне рівняння: 


у-90-9ч5 


де 9 - географічна широта місцезнаходження ЗСП. 
Запишемо кут 1 у вигляді синусоїдального рівняння: 


2 -2345 "Р. (4-81 
365 


(5) 


(6) 
де 4 - порядковий номер дня року (1 січня - 4- 1). 
Кут у; визначає зміну радіаційного випромінювання Сонця протягом доби 
нік 
Іят 
Їсд 3 (7) 


де Їм - ЧАС КОрекції місцевий, с. 
Час корекції /км Враховує фактичний місцевий час (літній або зимовий). Він 
вираховується за допомогою рівняння (|З: 


рас . 3 які 
їх 7356-1070. ках 2) (8) 


ДФ Ірк тах - ЧАС КОЛИ Сонце перебуває у зеніті, с. 

Кут у», дозволяє враховувати вплив обертання Землі навколо Сонця. Так, у місяці 
рівнодення (березень, вересень) даний кут становить 0 або л, та не впливає на час 
перебування Сонця над лінією горизонту (прийнято вважати, що такий час складає 12 
годин). В інших випадках час перебування Сонця над горизонтом збільшується (влітку) 
або зменшується (взимку). Кут у: розраховується наступним чином: 

9-т) 


Уза-т 
си (9) 
де т - порядковий номер місяця року. 
Коефіцієнт К, визначається за допомогою рівняння: 


М 1 
п - 
1ч-К, 5іп У, (10) 
Поправочний коефіцієнт ХК, рахується достатнім, якщо ЗСП знаходиться в 
північній півкулі, та від'ємним в південній півкулі. Він визначається відповідно до 


табл. 1. 
Таблиця 1. Залежність поправочного коефіцієнта від географічної широти 


Географічна широта 9, град | 60 50 40 30 20 10 0 
Коефіцієнт Ки 0,8 0,5 0,4 0,3 0,2 0,1 0 
а 0 ЗАГаб 
А шо В шо Р . 
Позначимо срИ срРИ Не Ки 5 у, на «54 тм У. 5 ук, 
срИ 


М зна аа зіп А Ки зіп Узки 
срИ 


Далі підставимо (2) в (1) та, розв'язавши диференційне рівняння теплопровідності 
для конвекційного і радіаційного теплообміну, отримаємо вираз: 
засн по В(Авіп(р' -- Сі) - С сов(р' -- Сп))-- Е(Авіп(н -- С1) - Н соз(и -- Ст)) 

Й АЧС , 1) 
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(і-ї) 
НЕО 


І 
де а 5 


Рівняння (11) не дає змоги визначити інерційну складову, що присутня при 
циклічній добовій зміні температури ЗСП. Тому необхідно звести рівняння до виду, 
який дозволяє виділити дану складову. Використаємо відоме тригонометричне 
рівняння Юніса. Таким чином, отримаємо вираз: 


зіп(р'-- Сі Ф)н В 


Р! 
ня цена; Я СТ У) Н 
УР"? 1 ре є ) ) (12) 


Р шо А В я "ТИВ УЖ віп Ук, Фр шо агс8іп А З 
8 С сб реє, 


2 2 


ТаТ, ї 0ЗАТ, 


Рівняння (12) є розв'язком диференційного рівняння теплопровідності (1), та дає 
змогу визначати циклічну добову зміну температури ЗСП. Враховуються інерційні 
властивості ЗСП, а саме: відмінність амплітуди температури поверхні ЗСП від 


-- -- АР? -1 
температури грунтової (асфальтованої) дороги визначається членом та 
/ 


УР он відставання температури поверхні ЗСП від температури навколишнього 
середовища визначається значенням кута м. 

Проведемо розрахунки, використовуючи математичну модель (12). Крім цього, 
використаємо дані, отримані за допомогою експериментальних досліджень, що 
проводились у період з 13 по 14 серпня 2013 року. На основі отриманих результатів 
побудуємо відповідну графічну залежність (рис. 3). 

Коротко опишемо основні параметри ЗСП. Для розрахунків брались наступні 
об'єкти місцевості: ЗСП без рослинності (рілля), ЗСП вкриті рослинністю, 
асфальтована дорога. У зв'язку з великими розмірами ЗСП, вважається, що вони 
розташовані на необмеженій площі, яка бере участь у теплообміні рівномірно по всій 
поверхні, а її розміри зводяться до ділянки 1Х1Х0,05 м. Основні теплофізичні і масо- 
габаритні дані ЗСП наведено в табл. 2. 

Присутня рослинність 


Т, . Вілля (експеримент). « ікелеримені) 
/ 


Присутня рослинність 
299 рисутня р 
297 


295 


ан 
золю. Асфальтована дорога 


293 


291 


Рілля 
289 Кр поки в о п и п и п п п п и и и п п я п п и п о р п о п и п и о п и в и в и о 
9.00 11:00 1300 1500 17:00 1990 2100 2300 100 300 500 700 990 


Рис. 3. Циклічна добова зміна різних об'єктів місцевості 
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Таблиця 2. Значення теплофізичних та масо-габаритних параметрів 


Фізичні з 9:00 13.08.2013 по 9:00 14.08.2013 
величини Рілля Присутня рослинність (| | Асфальтована дорога 
тек 295 295 295 
ЛТ,К 12 12 12 

і 9-3,6-103 9-3,6:103 9-3,6:103 
іга С 6-3,6-103 6-3,6-103 6-3,6-103 
НИ 14-3,6-103 14-3,6:103 14-3,6:103 
ізн б 12:3,6:103 12:3,6:103 12:3,6:103 
5, м? 1 1 1 

а0 0,95 0,9 0,6 

до Вт/м? 1366 1366 1366 

ЛЕ 0,36 0,36 0,36 

М. 0 0 0 

9, град 50,5 50,5 50,5 

а 214 214 214 

а, Вт/(м?:К) 15 19 15 

с, Дж/(кг- К) 1100 4218 1670 

р, кг/м3 1250 1000 1100 

У, м? 0,05 0,05 0,05 
ікнівоба 46800 46800 46800 

т 8 8 8 

Ко; 0,5 0,5 0,5 


Загальна тенденція по циклічній зміні температури ЗСП більш наближена до 
експериментальних даних (рис. 3), ніж у роботі (2). Характерною особливістю 
математичної моделі є те, що амплітуда зміни температури ЗСП знаходиться в межах 
289...302 К, що відповідає експериментальним даним. На рис. 3 присутній інерційний 
зсув температури ЗСП різних видів, що не є характерним для роботи |2|. Даний 
інерційний зсув зумовлений теплофізичними властивостями різних видів ЗСП. 

Висновки 

Таким чином, рівняння (12) являє собою математичну модель прогнозування 
циклічних добових змін температури ЗСП. При її розроблені, на відміну від 
математичної моделі в роботі |2|, були враховані додаткові фактори впливу на зміну 
температури ЗСП, а саме: географічна широта місця розташування ЗСП, пора та день 
року, місцевий час (літній, зимовий), хмарність, додаткові теплофізичні параметри 
(питома теплоємність, густина). Існує можливість визначати інерційні властивості ЗСП 
за допомогою кута м. Порівняльний аналіз математичної моделі (12), математичної 
моделі роботи |2| та експериментальних даних вказує на те, що математична модель 
(12) дозволяє отримувати більш наближені до експериментальних даних значення 
температури ЗСП. Отже, існує можливість прогнозувати час доби, коли використання 
тепловізорів раціональне (за рахунок максимальної різниці значень між температурою 
ЗСП та інших об'єктів місцевості). Це дозволить враховувати даний факт при 
плануванні часу використання БЛА. 

Перспективними напрямками подальших досліджень можна вважати: врахування 
впливу циклічних теплообмінних процесів у навколишньому середовищі на 
акумуляцію тепла ЗСП; формалізацію оцінювання ефективності застосування 
розробленої математичної моделі. 
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ВЕ5ОМКЕ 

М.У. КогобсеНіпз5Кіу, М.О. 5Іопоу, А.А. Магійу 

Мафетагіс піодеі ої Іогеса5біпо дайу сусіеє гетрегаїиге спапоее5з ої аргісийигаї 
Іапдз 

Тре обіесі ої Фе агіїсіе 15 Гогеса5іпє, дайу сусіє іептрегагиге ехсрапоез ої аргісиіита! 
Іапд5 їо іпсгеа5е гепабійсу ої плопіїогіпє гезиіїв муріїе ц5іпо цшптаппей аєгіа! уеПісіез 
(ОАУз) едшрредй міїв Фегта! ітаєег. 

Тре азвідптепі ої Ше агісіе 15 сгеайпе плаїбетайс тодеі ої Гогеса5іипя даїу сусіе 
їеппрегаїиге ехспапее5 ої аєгісиійига! Іапд5 сап Каке іпіо сопз8ідегайоп плоге Гасіог5 аї 
тайетайс подав! їп ракепі (21. 

Тре дмезійоп ої плпіїогіпє ої аєбгісийига! Іапд5 Ба5 а зресіа! ріасе їп їбе сигтепі 
деуеіортепі ої Фе аєгісийига! 5есіог. ТЬі8 аПом/8 іо геасі абоці Шеїг 51шайоп апа їо ргедїсі 
йииге србапеез їп уїе!48. Опе ої Ше сПпеарезі апа тобі еНеспує мауз ої мізца! плопіїогіпо ої 
аєгісиига! Іапд8 15 цзіпе ап ПАУЗ. ТПезе ОАУбз аге едиірреай умій ісіемівіоп 5игуеШапсе 
сатега8 орегайтя їп Ше у18їб1е апа іпітагед гапее ої еіесіготаєпейіс зресітит. ТРі8 іпітагей 
гапое раз 5еуега! адуапіарез раї аПоху птлоге еНесіїує плопіїогіпо ої аргісиита! Іападз. 

Тре едпайоп (12) 15 плафетайс тоде! ої Гогесазіте, даїйу сусіє кетрегаїиге спапее8 
ої аргісиига! Іапд5. Деуеіоріпє їбе плодеі, ууе сопз5ідегей ехіга Гасіог5 ої іпйиепсе оп 
гепарегакиге спапеез ої аєтісиНига! Іапд8, зисі и8: Іайиде ої аргісиига! Іапд5, зеазоп апа 
Чау ої Ше усаг, Іоса| їте, сіопдіпез5, ехіга іегпойзіса! рагатеїсг5 (зресійс іїБегтаї 
сарасіїу, деп5іїу). ТПеге 15 роз51біШсу їо деїегтіпе Шегпаї іпегіїа сфагасіегізіїс ої 
абгісиіига! Іапд8 ц5іпе м апеіе. Апаїузі8 ої плаїдетакйс подеі (12), Г2| апа ехрегітепі даїа 
15 5Помлтпє аб плаїбетайс тоде! (12) сап бтуе плоге шиїрійі гезиїв Гог ехрегітепі ага. 50, 
її епабіе8 їо Гогеса5і дау те убеп Шегта! ітаєег цзіпє 15 гайопа! (тахітит фіНегепсе ої 
їепттрегаїиге Бебмееп аєстісиита! Іапд5 апа обег обіесіє іп їБбе агеа). Тетрегаите Рогеса5і 
аПому8 Гог Бебег йте ріаппіпо уубеп цпзіпє ПАУЗ. 
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